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ABSTRACT 


Temperature, velocity and viscosity profiles for coupi»_d 
thermal and mechanical models of the upper mantle beneath 
continental shields and old ocean basins show that under 
the continents, both tectonic plates and the asthenosphere, 
are thicker than they are beneath the oceans. The minimum 
value of viscosity in the continental asthenosphere is 
about an order of magnitude larger than in the shear zone 
beneath oceans. The shear stress or drag underneath con- 
tinental plates is also approximately an order of magnitude 
larger than the drag on oceanic plates. Effects of shear 
heating may account for flattening of ocean floor topography 
and heat flux in old ocean basins. 


INTRODUCTION 


Plate tectonics has stressed the importance of the con~ 
cepts of lithosphere, the rigid surface layer of the Earth 
which is divided into a mosaic of moving plates, and of 
a stheno sphere, a softer layer under the plates facilitating 
the mechanical decoupling with the lower mantle. The con- 
cept of an asthenosphere was derived long ago from the study 
of isostasy. On the other hand, the physical structure of 
the upper mantle has been discussed either on the basis of 
seismic velocity distributions derived from observations 
(Press, 1972; Knopoff, 1972; Jordan, 1975), or on the basis 
of purely thermal models (McKenzie, 1967; Sclater and Franche- 
teau, 1970; Turcotte and Oxburgh, 1972) . Both of these latter 
approaches have led to the conclusion that the physical 
structure is different under oceans and under continents. 

It is, however, difficult to make quantitative comparisons 
between seismological data and theoretical therroal profiles, 
since this involves a precise description of the petrologi- 
cal state of the mantle. For oceanic plates a^one other 
geophysical data like topography have been used to test 
theoretical thermal models (Sorokhtin, 1973; Davis and Listen, 
1974; Parsons and Sclater, 1975; Oldenburg, 1976) . Convec- 
tive models have not focused on the distinction between 
oceanic and continental regions. 
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Here we present models of the upper mantle wherein 
both the thermal and flow velocity profiles are calculated 
on the basis of coupled thermal and mechanical equations 
satisfying simple geophysical boundary conditions (Froidevaux 
and Schubert, 1975; Schubert et al., 1976). This is illus- 
trated in Figure 1 which shows the oceanic model on the 
left, and the continental one on the right. The boundary 
conditions indicated in the figure, are that the temper- 
ature and horizontal velocity at the surface be T^ and u^, 
and that the horizontal velocity at great depth be zero. 

In the continental case the heat flows at the surface 

o 

and at great depth q^ are also imposed, whereas in the 
oceanic case,, the temperature at great depth T^ and the 
vanishing of the vertical velocity at the surface are 
specified.' The physical quantities to be calculated are 
also listed in the figure; they are temperature T, hori- 
zontal velocity u, shear stress t and viscosity, and, for 
the Oceanic case only, vertical velocity v. Except for 
T, these quantities vary with depth y. In the oceanic 
case they all vary with horizontal distance from the ridge 
X, i.e. with age x/u^. Given a petrologic model of the 
mantle , seismic velocity profiles can also be computed. 

All models are based on the law of deformation of 
olivine (Post, 1973; Kohlstedt et al., 1975) 
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where T is the absolute temperature, p is the pressure, 

e H ^ du/dy is the straxn rate, and R is the gas constant. 

The proportionality constant ^ and the activation energy 
* 

E have been measured in the laboratory, whereas the acti- 

•k 

vation volume V has to be estimated. The justification 
for a model based on olivine is that this mineral is the 
most abundant in the mantle and it deforms more readily 
than other minerals in most mantle-derived peridotite 
samples (Boullier and Nicolas, 1973). The coupling between 

velocity and temperature comes about through the T-dependence 

* 

of {!' and through the viscous heating term 2ie in the 
temperature equation. In the continental case the temper- 
ature equation also contains a radioactive heating term 
and a vertical conduction term. The continental mantle 
is overlain by a crust divided in 2 layers. Radioactive 
heat generation in the upper 8 km is a given quantity, 
whereas the volumetric heat production Q 2 for the 32 km 
thick lower crust has to be calculated to satisfy the boun- 
dary conditions. In the oceanic case the temperature equa- 
tion contains both horizontal and vertical advection of 
heat, vertical conduction of heat and viscous dissipation. 
Since horizontal advection of heat is a dominant effect, 
radioactive heating has been neglected in this presentation. 
The vertical velocity in the oceanic model requires consider- 
ation of another mechanical equation, the continuity equation. 
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These models imply a mechanical decoupling between the 
surface, moving at velocity u^, and deeper regions. They 
do not fix, a priori, the depth of the shear zone nor the 
magnitude of t. Since no return flow is assumed at shallow 
depth, the momentum equation simply states that t is inde- 
pendent of y; in the oceanic case it will however vary 
with X, the horizontal coordinate or age. 

The system of equations just described contains no 
time dependence. The models should thus apply in regions 
where steady state might exist. We assume this to be the 
case for oceans and for old shields (Polyak and Smirnov, 
1963) . The quasi-equality of heat flow for old continental 
shields and old ocean basins has been considered a paradox. 
Our models should tell us something about this problem. 

They should also lead to an interesting comparison of the 
lithcsphere-asthenosphere structures of both regions. Such 
comparisons will be presented in section 3, after the most 
striking characteristics of the models have been reviewed 
for continents in section 1 and for oceans in section 2. 
Further details and mathematical formulations are given in 
Froidevaux and Schubert (1975) and Schubert et al. (1976) . 
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1. TEMPERATURE AND VELOCITY UNDER OLD CONTINENTAL SHIELDS 

The values of the rheological parameters in (1) depend 
on the amount of volatiles# particularly water, present in 
the mantle (Carter and Av4 Lallemant, 1970; Stocker and Ashby, 
1973) . It is therefore useful to test the behavior of our 

* -k 

thermo-mechanical model when B^, E and V have values ap- 
propriate to wet or dry olivine. This is illustrated in Figure 

* 

2, where solutions are drawn for extreme values of V (11 and 

3 

30 cm /mole) for a continental system with mantle radioactivity 
but without heat flow from great depth. The hatched areas in 
Figure 2 indicate the range of solutions for intermediate val- 
ues of V . Other parameter values are listed in the iigure 
caption. The temperature rises monotonically to a final value 

4r 

depending mainly upon V . This dependence reflects the con- 
tribution of shear heating to the thermal state. The velocity 
drops from its surface value, 5 cm/yr, to zero in a shear 
zone or asthenosphere beneath the rigid plate or lithosphere. 
The depth of the shear zone is shallow or the lithosphere is 

thin when the pressure dependence in (1) is strong, i.e. when 
* 

V is large. A larger activation volume also yields a thinner 
asthenosphere. These characteristics are shown by the dotted 
bands representing the shear zones for the various solutions. 
The temperature and velocity profiles are relatively insensi- 
tive to whether the mantle is wet or dry. One remarkable 
property of all these solutions, is that the viscosity reaches 
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about the same minimum value of 10 poise in the shear zone, 
or asthenosphere . Thermo-mechanical coupling leads to self- 
adjustment of the temperature and of the depth of the astheno- 
sphere, so that the values of the rheological parameters do 
not influence the value of the viscosity minimum. 

What variables in our continental models have a major 
influence on the viscosity value in the asthenosphere? There 
are two: the amount of radioactive heat sources in the 

mantle and the plate velocity. The effects of radiogenic 
heating in the mantle are seen in Fiaure 3 which oresents 

continental solutions for u = 2 cm/vr with = 0 HGU and 

o m 

Q = 0.06 HGU (1 HGU = 10^^ cal/cm^ sec). The addition of 
m 

radioactivity raises the temperature at all depths, moves 
the asthenosphere upward and reduces the value of the minimum 
viscosity bv an order of magnitude. By comparing solutions 

* 3 

in figures 2 and 3 for V = 11 cm /mole and Q =0.06 HGU one 

m 

finds that an increase in plate velocity from 2 to 5 cm/yr 
lowers the viscosity by another factor of 4. The viscosity 
profile for V = 11 cm /mole, (^ = 0.06 HGU and = 5 cm/yr 
is shown as the dashed curve in Figure 3, The computed shear 
stress T, i.e. the drag under the continental plate, also de- 
creases with increases in mantle radioactivity and plate vel- 
ocity. For the cases in Figure 3 with u^ = 2 cm/yr, the t 

values are 114 bars for Q = 0 HGU and 27.2 bars for = 0.06 

m m 

HGU; in the latter case t drops to 17.1 bars when u^ is in- 
creased to 5 cm/yr. This decrease of t with u^ is not as 
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strong as that found by Froidevaux and Schubert (1975) in 
situations where viscous heating was the dominant effect in 
the mantle and tu^ was almost constant in the velocity range 
investigated. 

Comparisons with the Earth have been attempted and the 
results are shown in Figure 4. Since the activation volume 
is not known we imposed an additional constraint on the model# 
naunely that the geotherm is required to pass through the data 
points (shaded areas in Figure 4) derived from phase equilibria 
in kimberlite nodules from the South-African shield (Boyd# 

1973; MacGregor and Basu# 1974) . Only the unperturbed branch 
(Green and Gueguen# 1974) of this experimentally determined 
geotherm is used. Solutions for wet (Post, 1973) and dry 
(Kohlstedt and Goetze# 1974) olivine rheological parameters 
are obtained with q^ = 0.1 HFU# the imposed heat flow from 
the deep mantle. This value of q^ corresponds to an estimate 
of the adiabatic temperature gradient. The temperature solu- 
tions in Figure 4# which approximately satisfy the nodule 
data, are relatively insensitive to the value of the activa- 
tion volume- Thus the model geotherms alone do not provide 

* 

much of a constraint on the value of V . However# since the 
sheared nodules in kimberlites are thought to originate at 
depths of about 180 km or somewhat deeper (Boullier and 
Nicolas# 1973) # one can use the depth of the shear zone in 
the model solutions as an additional constraint on acceptable 
models and appropriate values of V . 
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* 3 

For wet olivine/ a V between 20 and 28 cm /mole would 
give models with acceptable temperatures and shear zones 
starting at depths between about 225 km and 180 km, respectively. 
For dry olivine, no solutions could be found which simultan- 
eously satisfied the constraints in both temperature and 
depth to the top of the shear zone (depth at which u/u^ =* 0.95). 
In the case of dry olivine V values of 16 and 20 cm /mole 
yield models with shear zones starting at depths of about 
305 km and 260 km, respectively. 

ThiS viscosity minima for both the wet and dry solutions 
have values of about 10^^ poise. The shear stresses for 
these solutions are abc , 80 bars, in agreement with values 
proposed by Goetze and Kohlstedt (1973) on the basis of dis- 
location studies in mantle derived peridotite nodules. The 
heat generation values Q 2 in the lower crust range from 0.75 
to 0.8 HGU, about 40% to 35% lower than averaged values for 
gabbro. In these solutions, the granitic upper crustal layer 
had Qj = 5 HGU and » 0.06 HGU was specified for the mantle 
down to 500 1cm depth. 

In conclusion, if a steady state regime is assumed to 
exist under old continental shields, our coupled thermo- 
mechanical model is capable of satisfactorily predicting 
the physical structure of the upper mantle. The finite 
thickness asthenosphere underlying the rigid lithosphere is 
a consequence of the strong T and p dependences of the mechan- 
ical constitutive equation (1) . Under old continents the 


I 


I 


I 


I 


4 


9 

♦ 

depth of the asthenosphere is determined mainly by V , whereas 
its ’softness' is essentially governed by and u^. Shear 
wave velocity profiles computed on the basis of solid state 
data for pure olivine and our model geotherms are flat. Phis 
agrees with seismological data and emphasizes that the seis- 
mological low velocity zone, which does not exist under old 
continental shields, is not equivalent to the shear zone or 
asthenosphere . 
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2. TEMPERATURE AND VELOCITY UNDER OLD OCEANIC BASINS 

For oceanic plate models, the two-dimensional flow field 
u and V and temperature profile depend not only upon depth but 
also upon distance from the ridge as pictured in Figure 1. 

The boundary conditions do not include the surface heat flow, 
since this quantity is not well known. Instead we have im- 
posed the temperature T^ at great depth. This quantity is 
not known either, so that several values between 1200®C and 1600°C have 
been used. This boundary condition has the advantage of 
Baking a comparison with purely thermal models simple. All 
models in this and the remaining parts of the paper are ba&ed 
upon a dry olivine rheology. 

The temperature profiles are dominated by heat advection, 
i.e. by lithospheric and asthenospheric cooling. This is 
particularly true for high T^ values and slow plates. For a 
rather cool asymptotic temperature = 1200**C and a high 
plate velocity u^ ■ 10 cm/yr. Figure 5 shows solutions for 2 
ages, 10 and 150 My. For 10 My the T profile is almost ident- 
ical to the purely thermal solution calculated for the same 
variable thermal conductivity (dotted curve) . The velocity 

prpfiles u (y) and v (y) are seen to depend on the value of 

* * 

V ; the shear zone is shallower when V is large, as was 

found in the continental models. In Figure 5 the hatched 
areas indicate again the range of solutions for intermediate 
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* 21 
values of V . The viscosity minimum is less than 10 poiso, 

for « 1600 C it drops below 10^ poise. The vertical 
velocity v increases from zero to an asymptotic value 'n the 
region where the horizontal velocity drops to zero, i.e. in 
the asthenosphere. The finite value of v at great depth 
provides a mass flow necessary for the accretion of the thick- 
ening lithosphere. For an age of 150 My Figure 5 shows that 
the temperature profile departs strongly from the purely ther- 
mal cooling solution (dotted curve) . This is due to shear 

heating, an effect which is stronger for larger values of V** 
The vis€ o- ity in the asthenosphere is somewhat higher be- 
neath oceans. At 150 My the upward velocities are an 
order of magnitude smaller than at 10 My. 

Shear heating generates marked departures from a simple 
/age dependence for derived quantities like lithosphere thick- 
ness, topography and heat flow. This is illustrated in Figure 
6 for the first two of these geophysical quantities. The 
lithosphere thickness, defined as the depth at which u falls 
to 95% of the surface value u^, increases with age, but not 
as fast as a /age low which would yield a straight line on 

this plot. As mentioned earlier the thickness is smaller 

* * 

for larger V values. For each value of V the upper and 

lower curves correspond to u^ values of 10 cm/yr and 2 cm/yr, 
respectively. Shear heating is stronger for larger u^ values; 
consequently the lithospheric plate is somewhat thinner. The 
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topography of the ocean floor represented in Figure 6 

by the water-depth also departs from a /age law indicated 

by the dashed straight line. The departure is very 

* 3 

pronounced for V = 30 cm /mole. The experimental 

points redrawn from Parsons and Sclater (1975) could be 

* 3 . 

fitted for a V value just above 11 cm /mole. This 

topography model reflects the gradual cooling of the lith- 
osphere and asthenosphere ; it was computed for a thermal 

-5 -1 

expauision coefficient a = 4x10 K . No data is avail- 
able for the lithospheric thickness as defined here. The 
depth of origin of sheared peridotite nodules might in 
the future give some constraints. For now we can only 
say that the depth to the seismic low velocity zone infer- 
red from surface and body wave data increases monotonically 
up to ages as old as 150 My (Leeds, 1975; Sipkin and 
Jordan, 1976). The ocean floor heat flow data, on the 
other hand are too controversial to warrant a detailed 

comnarison. Our models of q for T = 1200^0 also show 
^ o °° 

a certain amount of flattening, i.e., a departure from 

- 1/2 

an (age) ' proportionality for heat flow. Viscous 
dissipation provides additional heat to the old ocean floor. 
Stress values increase by factors of 2 to 5 for the cases 
shown above as age increases. They are less than 80 bars 
for solutions based on the rheology of dry olivine. For a 
wet olivine rheology, stresses and shear heating effects 
are reduced by factors of 2 to 3. They are also reduced for 
larger T^ values, = 1200°C being the value applicable to 
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Figures 5 and 6. 

In suiranary, for young ages thermo-mechanical coupling 
beneath oceans is negligible and the T profiles are close 
to the solutions for simple boundary layer cooling. Note 
however that although the isotherms for such purely thermal 
solutions are parabolic, they would not yield a parabolic 
lithosphere because of the pressure depend -ice of the vis- 
cosity implied by equation (1) , A second point to emphasize 
is that departures from /age laws at old ages are a result 
of distortion of the geotherms by shear heating. This ex- 
planation of known geophysical behavior (see topography) is 
more sat- 3 Tactory than that based on the ad hoc representa- 
tion of the lithosphere by a slab of constant thickness 
(Parsons and Sclater, 1975) . The stage is now set for a 
comparison of continents and oceans. 
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3. A DIRECT COMPARISON OF CONTINENTAL AND OCEANIC MODELS 
BASED ON A DRY OLIVINE RHEOLOGY. 

So far we have shown that the one-dimensional flow field 
under old continental shields, as well as the two-dimensional 
flow field under oceanic plates, and their associated temper- 
ature fields can be modelled on the basis of simple geophysical 
boundary conditions and physical properties corresponding to 
the thermal conductivity and rheology of olivine as measured 
In the laboratory. 

A comparison of the situations under old continents and 
under old oceanic basins can be undertaken. This requires 
an appropriate matching of boundary conditions. We shall pro- 
ceed in the following way. First we take T^ = O^C, = 1.05 

—6 2 

HFU (IHFU = 10 cal/cm sec) for continental shields and 
compute solutions corresponding to q^ = 0 HFU. This yields 
a Computed temperature T^ at great depth under old continents. 

We then compute oceanic solutions for the same T^ and T^ 
values, and of course for the scime surface velocity u^. All 
this is carried out under the assumption that mantle rheo- 
logical properties are identical to those of dry olivine 
and are the same under ocean basins and under continental 

shields. Since the activation volume V has not been measured 

* 3 

in the laboratory, we present solutions for ll^V <30 cm /mole 
In the oceanic models, where horizontal and vertical ad- 
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vection of heat dominates, we shall neglect the thin basaltic 
crust as well as mantle radioactivity. However the continental 
models not only include a radioactive crust with = 5 HGU 

and a lower crust with a Q 2 value which has to be computed, 
but also a radioactive upper mantle with = 0.06 HGU, a 
reasonable average value. We realize that the true value of 
is not well determined; for example Sorokhtin (1974) has 
argued that its value could be zero. In any case, although 
the specific comparison discussed here is based on a contin- 
ental model with nonzero Q we can use the results of Figure 3 

m 

to infer the effect of a smaller or zero value of beneath 
continents. Shear heating is, of course, included in both 
oceanic and continental models. 

Figure 7 gives continental and oceanic solutions for 
* 3 3 

u^ = 2 cm/yr. For V = 11 cm /mole and 30 cm /mole, the 
continental solutions have T = 1392 Q- = 1.04 HGU and 

00 ^ 

= 1665 °C, Q 2 = 0.59 HGU, respectively. The correspond- 
ing oceanic geotherms for an age of 150 My have surface heat 
flow values of 0.82and0.98 HGU, respectively. Both the 
surface heat flows in the oceanic cases and the lower con- 
tinental crustal radioactivity have reasonable values when 
compared with geophysical and geochemical data. 

What are the most striking points emerging from a com- 
parison of the solutions in Figure 7? Firstly, we see that 
the oceanic geotherm is everywhere higher than the corres- 
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ponding continental one ctnd thus has a much steeper grad- 
ient at shallow depths. This is basically due to the exist- 
ence of a thick radioactive crust for the continents. As 
a consequence the lithospheric thickness is larger under 
continental shields than under old oceanic basins. This is 
pictured by the position of the top boundary of the dotted 
bands representing the depth range of the asthenosphere 
(defined by 0.05 £ u/u^ £ 0.95). The asthenosphere beneath 
continents is thicker than that under oceans. The value of 
the viscosity minimum in the continental asthenosphere is 
essentially an order of magnitude larger than that in the 
oceanic shear zone. Similarly the drag under the continental 
plate is about an order of magnitude larger than that beneath 
the oceanic plate. The drag is given by the computed value 

of shear stress x which is 4.9 and 27.2 bars for the oceanic 

* 3 

and continental cases, respectively, when V = 11 cm /mole; 

* 3 

corresponding x values for V = 30 cm /mole are 10.7 and 
124 bars. These results compare favorably with those of 
Forsyth and Uyeda (1975) and Solomon et al. (1975) who 
suggest a larger drag under continental plates on the basis 
of analyses of plate motions. 

To discuss the relative importance of the various heat 
sources beneath the oceans and continents we have constructed 
Table 1 from the solutions of Figure 7. In order to make 
vertical heat advection into the oceanic lithosphere negligible 
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we now define the base of the oceanic lithosphere as the 

depth for which u = 0.995 u^. For consistency we apply the 

same definition of lithospheric base to the continental case, 

even though there is no advection in the continental models. 

Given this definition we can read the entries in Table 1 

for lithospheric thickness and sublithospheric temperature 

directly from Figure 7. The heat flux into the base of the lith 

osphere, q_, is computed as the product of thermal conductiv- 

ity and temperature gradient at the appropriate depth from 

the detailed solutions. The difference between surface 

heat flow and q^^ is either the contribution of lithosphere 

radioactivity to q^ in the continental case or the contribution 

of horizontal advection to q in the oceanic case. 

o 

Table 1 shows that a substantial fraction (1/3 to 2/3) 
of the surface heat flow in the old ocean floor comes from 
the cooling of the lithosphere as it thickens and moves away 
from the ridge. The remainder of the surface heat flow de- 
rives from a stheno spheric cooling, vertical advection, and 
viscous dissipation tu^. The source of the continental 
surface heat flux is radiogenic heating in the crust (55 to 
70%) and in the mantle (26%) , and viscous dissipation (19 
to 4%) . 

Sclater and Francheteau' s (1970) purely conductive 
model of the oceanic lithosphere does not include advection 
of heat. Our results show that this is still important even 
for an age of 150 My and must be included when the "equality" 
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of heat flow between continental shields and oceanic basins 
is discussed. One should notice finally that the sub-astheno 
spheric temperature is considerably higher (200-400°C) than 
the temperature just under the plates. In the literature 
these 2 quantities are often indiscriminately equated, just 
as the constant thickness slab introduced by McKenzie (1967) 
for mathematical convenience is often presented as the lith^ 
o sphere . 


CONCLUDING REMARKS 
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This paper has shown that it is possible to construct 
satisfactory models of the upper mantle beneath both 
continents and oceans. Under the oceans, the plates are 
found to be thinner, the asthenosphere softer and hence 
the drag smaller than under the continents. The comparison 
of continental and oceanic models with each other and with 
the real earth would be improved by including radioactivity 
under the oceans and by accounting for a nonzero heat flux 
from great depth. Further improvement would require addition- 
al data on possible lateral heterogeneities in the distribu- 
tion of radioactive heat sources and in the rheological 
properties of the mantle. 
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FIGURE CAPTIONS 


Figure J. . Oceanic (left) and continental (right) models. 

The boundary between lithosphere and astheno- 
sphere appears to be sharp for artistic reasons 
only; its depth is no'^ a specified parameter. 

For the oceanic case the temperature hori- 
zontal and vertical velocities u^ and v^ at 

o o 

the surface, as well as the temperature 
and the zero horizontal velocity u = o at 
great depth are given as boundary conditions. 

For the continental case the boundary conditions 
are the temperature T^, the velocity u^ and the 
heat flow at the surface, as well as the 

heat flow and zero velocity u = o at great 

depth. The continental mantle has a volumetric 
radiogenic heat-source Q^. It is overlain by 
a crust with thermal conductivity Ic . The 
crust has an 8 Ian upper layer with radiogenic 
heat sources and a 32 )cm lower layer with 
unspecified radiogenic heat sources Q 2 . The 
quantities listed on the right hand side of 
both the oceanic and continental diagrams are 
those to be computed. 
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Figure 2. Temperature and velocity profiles for a wet and 

a dry olivine continental mantle containing 

radioactive heat sources. The parameter values 

-15 3 5 

for the wet case are B., = 5.41x10 cm sec K, 

^ 3 

* 

E =95 kilocal/mole (Post, 1973) , and for 

the dry case B, = 6.45x10 cm^sec^K, E = 125 
3 . . 3 

gm 

kilocal/mole (Kohlstedt and Goetze, 1974) . The 

* 3 

activation volume V was 11 or 30 cm /mole and 
the corresponding solutions are labelled accord- 
ingly, The hatched area indicates the range of 

■k 

solutions for intermediate values of V . In all 
of these figures, the mantle thermal conductivity 

model k (T,z) of Schatz and Simmons (1972) is 

m 

used. The continental crust has a top layer of 
8 km with = 5 HGU and a bottom layer of 32 km 
with the computed value quoted :'n the text. 

The mantle radioactivity content is 

specified down to 500 km depth. Boundary con- 
ditions are T = 0 °C, q =1.05 HFU, u = 5 cm/yr, 
o o o 

= 0 HFU, The computed shear stress t is 9.3 bars 

(V =11 cm^/mole) or 50.9 bars (V = 30 cm^/mole) 

* 3 

for the wet case, and 17.1 bars (V = 11 cro /mole) 

* 3 

or 63.2 bars (V = 30 cm /mole) for the dry case. 

The dotted bands indicate the depth range of the 
SSthenosphere for a particular value of V in cm ’'/mole. 


1 


I 


2: 

The asthenosphere is defined as the depth range 

between u = 0.95 u and u = 0.05 u . 

o o 

Figure 3. Temperature, velocity, and viscosity profiles for 

a dry olivine, continental mantle. Q =0.0 HGU or 

m 

0.05 HGU, V =11 cm^ mole”^, u^ = 2 cm yr”^, and 
q =0.0 HFU. For the models with mantle radio- 
activity, the heat sources extend down to 500 km 
depth. The computed shear stresses are 27.2 bars 

(Q^ =0.06 HGU) and 114 bars (Q^ = 0.0 HGU). The 

* 

dashed curve is the viscosity profile for V = 

3 

11 cm /mole, Q =0.06 HGU and u =5 cm/yr (x = 
m o f j ' 

17 . 1 bars) . 

Figure 4. 'Realistic* models of the upper mantle beneath 

continental shields. The shaded region delineates 
the approximate range of the pyroxene geotherms 
(MacGregor and Basu, 1974). Computed models of 
the thermal and mechanical structures are shown 
for a range of activation volumes for wet and 
dry olivine. Geotherms satisfying the paleotemper- 
atures are relatively insensitive to the value 
of activation volume. An additional mechanical 
constraint based on the interpretation of sheared 
nodules (Boullier and Nicolas, 1973) is imposed 
on the depth of the top of the shear zone. For 
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wet olivine and q =0.1 HFU, the activation 

^00 

volumes for the solutions which best satisfy 

both the thermal and mechanical constraints 

3 -1 

lie between 20 and 28 cm mole . For dry 
olivine acceptable solutions could not be 
found. The computed shear stresses are 88.8 

bars (V* = 28 cm^/mole) , 71.9 bars (V =20 

3 * 3 

cm /mole) , 62.5 bars (V = 6 cm /mole) for 

* 3 

wet olivine, and 88.1 bars (V = 20 cm /mole) 

* 3 

and 62.4 bars (V = 4 cm /mole) for dry olivine. 

Figure 5. Temperature, viscosity, horizontal and vertical 
velocity fields of an oceanic lithosphere and 
asthenosphere model for a dry olivine rheology. 

T = 1200°C and u = 10 cm yr“^. Solutions 

CO O 

* 

are shown for two extreme values of V , 11 and 
3 -1 

30 cm mole , and for two ages, 10 and 150 My. 
They are indicated in the figure by the numbers 
adjacent to the curves. For an age of 10 My, 

temperature profiles are insensitive, at the 

* 

accuracy of our graph, to the value of V . 

This is not so for the velocity curves. The 

shaded area between two such curves represents 

the range of solutions for intermediate values 
* 

of V . For an age of 150 My the temperature 

★ 

profiles depend on V . The upper curve is for 
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* T — 1 

V s= 30 cm mole and the lower one is for 

* 3-1 , . 

V = 11 cm mole . The dotted curves denote 

the temperature profiles for simple boundary 
layer cooling with variable thermal conduct- 
ivity. The conduction solution for 10 My is 

* 3 -1 * 

similar to the V = 11 cm mole and V =30 

3 -1 

cm mole curves. The stresses developed in 

* 3 

these oceanic models are 11.8 bars (V = 11 cm /mole), 
34.0 bars (V* = 30 cm^/mole) for 10 My and 25.8 
bars (V* = 11 cm^/mole) , 75.5 bars (V* = 30 cmV^ole) 
for 150 My. 

Figure 6. Bathymetry of the ocean floor and lithospheric 

thickness (depth at which u/u^ = 0.95) as a 

function of (age) for a dry olivine rheology 

with T = 1200°C, and V =11 and 30 cm^ mole 

For the depth of water plot, u^ = 10 cm yr”^ 

-5-1 

and a = 4x10 K . The shaded region represents 

the bathymetry solutions for intermediate values 

* 

of V . The dashed lines are the ocean floor 
depths for simple boundary layer cooling. The 
data points are taken from Parsons and Sclater 
(1975) . For the thickness of the lithosphere 

curves, the value of distinguishes the solid 

* 

curves with the same value of V ; the thicker 
lithosphere is for u^ = 2 cm yr 


9 


the thinner 


30 


one is for u = 10 cm yr 
o 

Figure 7. A comparison of the thermal and mechanical 

structures between old oceanic basins (150 My) 

and continental shields. Dry olivine rheology, 

u = 2 cm yr = 0.0 HFU and equal asymptotic 

temperatures, T^, are required for both tectonic 

models. = 0.06 HGU is specified for the 

continental mantles down to 500 km depth. The 

* 

effect of an increase in V is greater on con- 
tinental models, which decouple at a greater 
depth than the oceanic ones. The viscosity 
minima for the oceanic cases are an order of 
magnitude lower than the continental values. 

The stresses in the continental cases (124.0 

^ # 3-1 * 

bars, V = 30 cm mole ; 27.2 bars, V =11 

3 -1 

cm mole ) are an order of magnitude greater 
than those found for the oceanic solutions 
(10.7 bars, V =30 cm^ mole 4.9 bars, 

V* = 11 cm^ mole'^) . 






















